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We investigate the twist-3 light-cone distribution amplitudes (LCD As) ol the scalar mesons ao, 
Kq and /o within the QCD sum rules. Lhe QCD sum rules are improved by a consistent treatment 
of the sizable s-quark mass effects within the framework of the background field approach. Adopting 
the valence quark component (9102) as the dominant structure of the scalar mesons, our estimation 
for their masses are close to the measured ao(1450), /To (1430) and /o(1710). From the sum rules, 



we obtain the first two non-zero moments of the twist-3 LCDAs 
and (6^*0) = 0.203 (0.093); those of the twist-3 LCDAs < 

},) = 0.018 (0.207); and those of the twist-3 LCDAs 



bag 



is. a 



(f K2) 



is,/? 



1(2) „) 



(£s,a ) 



0.369 (0.245) 
0.004 (0.355) and 

0.335 (0.212) and 



(£ct /q) = 0.196 (0.088), respectively. As an application of those twist-3 LCDAs, we study the 
B — > S transition form factors by introducing proper chiral currents into the correlator, which is 
constructed such that the twist-3 LCDAs give dominant contribution and the twist-2 LCDAs make 
negligible contribution. Our results of the B — > S transition form factors at the large recoil region 
q 2 ~ are consistent with those obtained in the literature, which inversely shows the present 
twist-3 LCDAs are acceptable. 

PACS numbers: 14.40.-n, 12.38.Aw, 11.55.Hx 



I. INTRODUCTION 

Even though lots of works have been done in the liter- 
ature, the properties of the light scalar mesons are still in 
ambiguity. In order to get an accurate theoretical predic- 
tion on the properties of the scalar mesons and on their 
applications to high energy processes, it is very impor- 
tant to provide a good interpretation of their complicated 
nonperturbative nature. 

Among the scalar mesons' non-perturbative sources, 
one of the most important thing is their light-cone dis- 
tribution amplitudes (LCDAs). At the present, some pi- 
oneering works for both the twist-2 and twist-3 LCDAs 
of the scalar meson have been done within the QCD sum 
rules, c.f. Refs.P, Hf. According to our experience on 
the light pseudoscalar twist-3 LCDAs, e.g. the pion and 
kaon electromagnetic form factors [!, |j| and the B — >• tt, 
K transition form factors @, Q, a well-behaved pseu- 
doscalar twist-3 LCDAs in the end-point region can give 
the conventional power suppressed contributions to the 
high energy processes in comparison to those of twist- 
2 LCDAs. It is interesting to know whether the scalar 
twist-3 LCDAs also possess such good feature. Moreover, 
a better understanding of the twist-3 LCDAs is crucial 
for a reliable estimation. The forthcoming more precise 
data, e.g. at the large hadronic colliders and the pro- 
gramming super B factories, also requires a more accu- 
rate theoretical estimation for the twist-3 contributions. 



' wuxg@cqu.edu.cn 



At present, we will investigate the twist-3 LCDAs of 
the scalar mesons by incorporating such quark mass ef- 
fects properly so as to achieve a more accurate theoretical 
prediction. Within the QCD sum rules, it has been found 
that the contributions from the quark mass terms (espe- 
cially those of the s-quark) will be comparable to that of 
the dimension-six operators or even the dimension-four 
operators, and it can even change the relative impor- 
tance of the operator expansion series counted by the 
naive power counting rules, c.f. Refs.[7Hl4j for study- 
ing the S'J7(3)-breaking effects of the kaon LCDAs and 
Refs.[lMl for the cases of the vector twist-3 LCDAs. 

For the purpose, we will calculate the Gegenbauer mo- 
ments of the twist-3 LCDAs for the scalars ao, Kq and /o 
within the QCD sum rules together with the QCD back- 
ground field approach. Basic assumption of QCD sum 
rules is the introducing of nonvanishing vacuum conden- 
sates such as the dimension-three quark condensate (gch , 
the dimension-four gluon condensate (G 2 ), and etc. [19j . 
The QCD background field approach provides a system- 
atic description for those vacuum condensates from the 
viewpoint of field theory (20l - [29j . It assumes that the 
quark and gluon fields are composed of the background 
fields and the quantum fluctuations around them. The 
vacuum expectation values of those background fields de- 
scribe the nonperturbative effects, while the quantum 
fluctuations represent the calculable perturbative effects. 
To take the QCD background field theory as the starting 
point for the QCD sum rules, it not only shows a distinct 
physical picture but also greatly simplifies the calcula- 
tion due to its capability of adopting different gauge con- 
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ditions for quantum fluctuations and background fields 
respectively. Because of the influence from background 
fields, the quark and gluon propagators shall include non- 
perturbative component inevitably, and the quark mass 
effect can be introduced in a consistent way. 

Moreover, the B — > S transition form factor within 
the light-cone sum rules (LCSR) provides a good plat- 
form for checking the properties of the scalar LCDAs. 
In the LCSR approach, a two-point correlation function 
is introduced and expanded near the light cone x 2 = 0, 
whose matrix elements are parameterized as LCDAs of 
increasing twists. By using the conventional currents in 
the correlator, the form factors will always contain the 
twist- 2 and twist-3 terms simultaneously, both of which 
play important roles for the final LCSRs f]. Because 
both the twist-2 and twist-3 LCDAs have their own un- 
certainties, the entanglement of them make the estima- 
tion under large uncertainty. Thus, for the sake of a 
better accuracy, it is helpful to choose proper chiral cur- 
rents in the correlators such that either the twist-2 or 
the twist-3 terms make no contribution to the LCSRs. 
In Ref. [30( . a chiral current has been suggested to make 
the twist-3 terms give zero contribution. At present, we 
shall introduce another type of chiral current such that 
the twist-2 terms make no contributions. Furthermore, 
we will deal with the semileptonic decays B — > Slvi and 
B — > Sll, which when in comparing with the forthcoming 
data shall be helpful for acquiring valuable information 
on the twist-3 LCDAs of the scalar particles. 

The remaining parts of the paper are organized as fol- 



lows. In Sec. II, we present the calculation technology for 
deriving the sum rules for the twist-3 LCDA moments of 
the scalar mesons ao, Kq and fo. Then, we present the 
formulas for the B — » S transition form factors. Numer- 
ical results and discussions for the scalar mesons' twist-3 
LCDAs and the B — > S transition form factors are given 
in Sec. III. The final section is reserved for a summary. 
In the appendix, we put the some more subtle points in 
deriving the sum rules for the twist-3 LCDA moments. 

II. CALCULATION TECHNOLOGY 
A. Twist-3 LCDAs of the Scalar Meson 




FIG. 1. Feynman diagrams for calculating the moments of 
the scalar meson's twist-3 LCDAs <fi s s and <j>g, where the back- 
ground gluon fields are included in the Fermion propagators 
implicitly and the background quark fields are depicted as 

crosses. The left big dot stands for the vertex ^iz ■ and 

(Jfiv (iz ■ D( n+1 '^ for cj> s s and (f> a s respectively. 

We adopt the suggestion of the valence quark contents 
dominant for the scalar mesons ao, Kq and /o, which are 
(du), (us) and (ss) respectively. The twist-3 LCDAs of 
the scalar meson are defined as Fy] 



(S(p)\q 2 (y)qi(x)\0)=m s fs due^+^^u) , (I) 

Jo 

(S(p)\to(y)v^q 1 (x)\Q)=-m s (p l >z u -p v z ll )fs [ due^+^^P-, (2) 



where z = y—x, ms andp are mass and momentum of the 
scalar meson, fs is the decay constant of the scalar meson 
defined by (S(p)\q~2<li |0) = msfs, u is the momentum 
fraction carried by q 2 quark, and u = l — u. The moments 
of twist-3 LCDAs are defined as 



(£) = / du(2u-l) w ^(«), 
Jo 

(o= [ du(2 U -irr s (u), 

Jo 



(3) 
(4) 



which satisfy 

(0|gi(0)(w- D) n q 2 (0)\S(p)) = m s f s (p ■ z) n {£) 

(mm^-DT+^^msip)) = 



i 



m s fs{PnZv-PuZn){p- z) n (C). (5) 
Here the zeroth moments have been normalized to one, 

(O = (&) = i. 

Fig.([T]) shows the Feynman diagrams for calculating 
the LCDA moments, where the background gluon fields 
are included in the fermion propagators implicitly and 
the background quark fields are depicted as crosses. 

To study the properties of these LCDAs, one can in- 
troduce the following two correlation functions, 
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i I dW*-<0|T{gi(x)(i*. Drq 2 (x) ) q 2 (0)q 1 (0)}\0) = -(z-q) n l^(q 2 ), 
i I d 4 xe i ^(0\T{q 1 (x)a^(iz- D) n + 1 q 2 (x), g 2 (0)<?i(0)}|0) = i{q^z v - q v z„){z ■ q) n I^°\q 2 ). (6) 



r 



Following the standard QCD LCSR within the back- we should deal with the sizable s-quark mass effects in 

ground field theory [20M22I . l24l [H, [2{| , we can derive the a more consistent way, which might cause sizable effects 

sum rules for the moments of <p s s and <f>g. For conve- comparable to those of the dimension-six condensates, 

nience, we present the detailed processes in Appendix A. For the purpose, we adopt the following propagators to 

During the calculation, as has been argued by Ref.[|j], do our calculation Q 

I 



S F (x,0)=i l -^e 



d 4 q - iq . x ( _ m+ j i J {j-m)^ 



m 2 — q 2 



(m 2 — q 2 ) 2 



bouL 



(m 2 — q 2 ) 3 ' (m 2 — q 2 ) 2 



(7) 



for the quark propagator and 
d 4 q 



iq-. 



■x J cab 



h«»\p = f[G^;p(0) + G w; ,(0)]. Here G V(t {x) = 
g s T a G1 jl (x) , the gauge invariant function G„ M;p (0) = 

g s T a Dp b G'l, fi (x)\ x= o, and the symbol ••• stands for the 
(8) irrelevant terms that lead to higher-order operators over 
than dimension-six. 

The sum rule for the even moments of (f> s s up to 
for the gluon propagator, where &oi^ = §Gv M (0) and dimension-six condensates is 



2,1 



-(2n + 3)x(l-a:) + 



m\m 2 — fn\ 2 
M 2 



M 4 e MMi-* )( fa, 



-in 







M 4 e- S °/ M2 dx 



2m\m 2 — m 2 2 



M 2 x(l - x) 



+ (a s G 2 ) f (2x-l) 2n ±- _(2n + l) + 
Jo 8?r . 



2m 2 + (2n + l)rni 2n(4n + \)m\ + 12nm\m 2 + 3mim 2 m\m 2 {2nmi + m 2 ) mfmf 



+ {g s q 1 aTGqi 



2 6M 2 2M 4 

9(2n — l)m 2 + n(16n — 5)mi ((8n — 3)mi + 3m 2 )m 2 m\m\ 



6M 6 



18M 2 



+g 2 s(qm} 2 



-8n 2 - 14n +12 m 



81M 2 



+.9s(<7i<7i)(<72<72); 



9 m\ — m\ 
+ [qi ^ ?2,mi ^ rn 2 ] >, 



27M 4 



-m\/M 2 



12M 4 



e- m V M2 +g 2 s (q iqi )(q 2 q 2 )~ 



9M 6 

2 /nf 2 



-ml/M 2 



' 2 '"'" : ' +\{ e -™W -I) 
m 2 



M 2 



-m%/M 2 ^ ■ « I e 



2m\m, 2 



-m 2 1 /M 2 _ 1 e -rr4/M 2 _ 1 



(9) 



where m 2 2 — m\x + m|(l — a;), mi and m 2 denote the masses of the qi and q 2 quarks, M is the Borel parameter. 
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The sum rule for the even moments of <^f s up to dimension-six condensates is 



imi/j e - m|/M2 <e> 



= I dx(2x - l) 2 "M 4 x(l - x)e ^»d-») 



i/^- 1 » 2 "f< 1 - I »( 1 + ^) + ^} M 



4 g -S <T /M 2 



(a s G 2 ) [ 
Jo 

(QiQi)- 
81 



24tt 
3mi + 



2mim2 
M 2 x(l - x) 



3 f/^3\ m l m 2 



/' 

JO 



, (2a; -l) 2 " 

dx \ ' — re m^ci-z) 



'2M 4 i(l -x) 



(4n + l)m\ m\m 2 



M 2 



+ 



M 4 



e -mlfM- + {gs q 1(TTGqi) 



24tt 2 



(16n + l)mi + 6m 2 m\m\ 



-An + 5 2m| 



M 2 



+ 



M 4 



36M 2 



+ 



9M 4 



-m^/M 2 



(10) 



The sum rule for the odd moments of </>| up to dimension-six condensates is 



\2n+l 



-2(ra + 2)x(l -x) + 



-2(n + 2)x(l-x)(l + -^] + 



l 2 L {2x - i] 
I! i)2n+1 

+ (a s G 2 ) / (2x-l) 2 " +1 ^ [_2(n+l) + 
+ (.9 S 3 /G 3 > 



171x1712 — 7Tl\ 2 
M 2 



M 4 e M-^ii-*)dx 
(2n + 3)m 2 2 + m\m 2 



2m\m 2 — m\ 2 



M 2 x(l - x) 



M 2 

2 

_ 12 

e M2 *(i-*> <£r 



M 4 e -5 3 /M 2 da , 



, 3 2n(2n+l) Z" 1 



(2a- l) 2 ™" 1 -^ dx 



+| - ('/...'/I 



96tt 2 

2m 2 + 2(n + l)m\ (2n + I) {An + 3)mf + 6(2n + l)m\m 2 + 2>m\m\ m\m\{{2n + l)m 1 + m 2 ) 



6M 2 



2M 4 



+ 



6M 6 
+5s(9i9i) 2 



e -ml/M~ + {gs g iaTGqi) 



36nm 2 + (2n + l)(16n + 3)mi ((8n + l)mi + 3m 2 )m 2 m\m 2 



36M 2 



+ 



+ 



9M 6 



2/M 2 



-2{2n- l) 2 - 14n + 5 m 2 . 



81M 2 



27M 4 



12M 4 



M 2 



+ ^2(e- m2/M -1) 



2 1 

+4?ra s (gi gi ) (q 2 q 2 ) - — 2 

- [qi «• g 2 , mi m 2 ] |. 



( e -</ M2 -e- m *l M2 ) + 2m im2 



-m 2 /M 2 _ 1 e -m 2 2 /M 2 _ ^ 



rrii 



(11) 



The sum rule for the odd moments of <j)g up to dimension-six condensates is 



-\m 2 s P s z- ml/M \£ n+1 ) 

J dx(2x - l) 2rl+1 M 4 x(l - x)e - — J dx{2x - 1) 



2n+l 



5 



1 



,2n+l 



9 2 {qlqlf 

81 



6 



24tt 

3mi 4 
-4n + 1 

+ 



(An + 3)mf m\m 2x 



M 2 
2m\ 



+ 



2mim2 
M 2 x(l - i) 

2 
2 



, (2a;-l) 2 "+ 1 

dx- ~. — ~, — re m2x ( i-* ) 



2M 4 x{l 



M 4 



M 2 



M 4 



l/M 2 



e ~m 2 /M + {gs q iaTGqi) 

[<?1 «?2,mi-H m 2 ] 



(16n + 9)toi + 6to2 m\m 



36M 2 



9M 4 



(12) 



B. B — > S Transition Form Factors within the 
Light-Cone Sum Rules 

The B — » S transition form factors f±{q 2 ) and 
/t(<? 2 ) are denned through the hadronic matrix elements 
(S(p)\q 2 7 l il5b\B {s} (p + q)) and (S(p)\q 2 a^-f 5 q"b\B {s) (p + 
q)), which are 



(S{p)\q 2l ^b\B {s) (p + q))=-2ip^ + { q 2 )-i[f + (q 2 ) + f-{q 2 )]q^ 



(S{p)\q 2 <J^ l5 q"b\B {s) {p + q)) = [2p„q 2 - 2? M (p- q)] 



-Mq 2 ) 
m B( ^ + m s ' 



(13) 
(14) 



These form factors are key factors for studying the 
semileptonic decays B/ s \ — > SWi and B/ s ) — > Sll. 

By introducing proper chiral correlators, we obtain the 
light-cone sum rules for those form factors that depend 

I 



I 

only on the twist-3 LCDAs of the scalar mesons, and then 
our twist-3 LCDAs derived in the last subsection apply. 
More explicitly, we suggest to calculate the following cor- 
relators : 



n M (p,g) = i / ^{SWlTfe^fl- 7s )b(x),b(0)i(l + 76)9i (0)}|0>, 
U»(p,q) = i I d i xe^^S{p)\T{q 2 {x)o^{l+ lb )q v b{x)M^+lMm^)- 



(15) 
(16) 



where qi,q 2 denotes the light quark field. Following the 
standard procedures to deal with the correlators which 



J 



are similar to that of the B — > pseudoscalar transition 
form factors, c.f. Refs. [3ll - l35j |. we can obtain the light- 
cone sum rules for the form factors /±,t(<7 2 ): 



/ + (9 2 ) 



f + (q 2 ) + f-(q 2 ) 



fAq 2 ) 



mb + m qi 



i B Jb 

msfs 
~6uM 2 



exp 



<t>s(u)(m 2 b - u 2 p- 



m b + m qi 
2 — o t-^— exp 



01 f B "l 



P s {u){ml+uV-<?) -^We X p(-A) , 



m s fs 
' 6u 2 M 2 



(m b + m qi )(m Bqi + m s ) 




uup — uq 



uM 2 



(j> s (u )exp 



msfs(u<f> s (u) + -<f>%(u)) 



so \ m 
M 2 ) m 



2 2 



+ uup — uq 
uM 2 



u\p 2 



ms/s(f s («) + ^^sW) 



(17) 



(18) 



du 



■ exp 



uup — uq 



ii 



M 2 



m b msfs , a , , 

x —<bo(u) 

3uM 2 Vsy ' 



ml+<p2 _ q2 )- (19) 



with 



u 



s - q 2 - p 2 ) 2 + %> 2 (m 2 - <? 2 ) - (s - q 2 - £> 2 ) 



III. 



2p 2 



NUMERICAL RESULTS AND 
DISCUSSIONS 



We adopt the follo wing input parameters to do our 
numerical analysis [37], l38| : 

m u = (2.5 ± 0.8) MeV, m d = (5.0 ± 0.8) MeV, (20) 
m s = (101 ± 12) MeV, a s (m c = 1.27 GeV) = 0.39,21) 



and 



(a s G 2 ) 
(9 3 sfG 3 ) 
(uu) 

(ss) 



(7.5 ±2) x 10~ 2 GeV 4 , 
(8.2 ±1) x (a s G 2 ), 
(dd) = -(0.254 ± 0.015) 3 GeV 3 , 
(0.74 ±0.03) x (uu}, 
(g s uaTGu) = (g s daTGd) = mg(w), 
ml = (0.80 ±0.02) GeV 2 , 
g 2 (uu) 2 = g 2 (dd) 2 = 2.693 x 10~ 3 GeV 6 , 
g 2 s (ss) 2 = 0.74 xg 2 (uu) 2 . 



(22) 



These condensates are given at the energy scale [i = 2 
GeV, which when necessary, will be run to the required 
scales by using the evolution equations [HI, [4(| ■ 



TABLE I. Masses for the scalar mesons. The threshold 
parameters are taken as Ss° a = (6-0 ± 0.5)GeV 2 , 5 Sj ° = 
(5.4 ± 0.3)GeV 2 and S f s % = (6.5 ± 0.3)GeV 2 . 



mesons 


M 2 (GeV 2 ) 


m(GeV) 


do 


1.037 ~ 1.666 


1.312 ~ 1.571 




0.991 ~ 1.501 


1.328 ~ 1.514 


fa 


1.284 ~ 1.803 


1.563 ~ 1.706 




1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.E 



M 2 (GeV 2 ) 



FIG. 2. Masses of scalar mesons do, Kq and /o versus M 2 
within their Borel windows. 



A. Masses and decay constants for the scalar 
mesons 

Setting n — in the sum rules (0 ITUl) for the moments 
of (f> s s and (j)^, we can further derive the sum rules for 
the masses of the scalar mesons. The mass sum rules 
are derived by doing the logarithm of these equations 
and applying the differential operation M A d/dM 2 over 
them. It is found that the results from both sum rules 
(|9l ITDl are consistent with each other, so we adopt the 
results from sum rule (J9j> for a detailed discussion. 

Usually, the threshold parameter S s (or S a ) is taken 
to be around the squared mass of the scalar's first ex- 
cited state. We adopt = (6.0 ± 0.5)GeV 2 for a Q, 

Sfj = (5.4 ± 0.3)GeV 2 and S{% = (6.5 ± 0.3)GeV 2 for 
Kq and /o 0] ■ According to the SVZ sum rule, the Borel 
window for the parameter M is determined by the re- 
quirement that the dimension-six condensate contribu- 
tion (SIX) does not exceed 10% and the continuum con- 
tribution (CON) is not too large, i.e. less than 30% of 
the total dispersive integration. In Table Q] and Fig.©, 
the Borel window is obtained by setting SIX < 1% and 



CON < 30%. The masses for the scalar mesons are pre- 
sented in TableUand their values versus M 2 are presented 
in Fig.©, whose central values are 

m ao = 1442MeV, m K * = 1421MeV, m fo = 1634MeV, 

These values are close to the physical states ao(1450), 
#£(1430) and / (17_10) [37|. This shows that the valence 
quark constituent du, us and ss are viable choices for 
studying the properties of a , Kq and f . 

Taking the scalar masses as inputs, we can further cal- 
culate the decay constants of the scalar mesons. For the 



TABLE II. Decay constants for the scalar mesons at the en- 
ergy scale fi = 1 GeV. 



mesons 


M 2 (GeV 2 ) 


/(GeV) 




1.6 ~ 2.5 


0.374 ~ 0.377 




1.5 ~ 2.2 


0.357 ~ 0.359 


fo 


1.8 ~ 2.2 


0.374 ~ 0.378 
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0.4- 



0.3- 



0.0-| 1 1 1 1 , 1 , 1 , 1 , 

1.4 1.6 1.8 2.0 2.2 2.4 2,6 

M 2 (Ge\f) 



FIG. 3. Decay constants of scalar mesons ao, Kq and fo 
versus M 2 within their Borel windows. 



TABLE IV. Moments from the tensor sum rules. 



mesons 


(&> 


M 2 (GeV 2 ) 




(&«„> = (0-203 ~ 0.204) 


2.0 - 2.5 


(Cla n ) = (0.092 ~ 0.094) 


2.0 - 2.5 




(£,*•> = (0-99- 2.54) xlO- 2 


1.8 - 7.8 


(£,*•) = (0.206 - 0.208) 


1.8 - 2.3 


fo 


(£*,) = (0.192 - 0.199) 


1.5 - 2.2 


(£, fn ) = (0.085 - 0.091) 


1.8 - 2.2 



purpose, the Borel window for each meson is redeter- 
mined following the same criteria as above, whose val- 
ues are collected in Table [TTJ The decay constants at the 
energy scale 1 GeV for the scalar mesons are presented 
in Table HI] and their values versus M 2 are presented in 
Fig-©; whose central values are 

f ao = 375MeV, f K . = 358MeV, f fo = 376MeV. 

In the following subsections, all the values of decay con- 
stants and moments are given at the scale 1 GeV unless 
explicitly pointed out. 

B. Moments for the scalar mesons 



FIG. 4. Moments of scalar meson ao versus M 2 within their 
Borel windows. 



and fo meson, and the first and second moments for Kq 
meson. Similar to the way to get masses and decay con- 
stants, we first find the stable Borel window for the sum 
rules of each moments. The choice of Borel parameters 
for each meson, together with the moments for the scalar 
mesons, are collected in Tables [TTT1 and HVl In the Tables 
IIIII and IIV1 the Borel windows are determined by setting 
SIX < 5% and CON < 30%. The moments for the scalar 
mesons versus M 2 are presented in Figs.(|H [6]). It is 
noted that the moments within their Borel windows are 
very stable and their central values are 



TABLE III. Moments from the scalar density sum rules. 



mesons 


<£.) 


M 2 (GeV 2 ) 


ao 


= (0.367 - 0.371) 


2.0 - 2.5 


(Clan) = (0.24 - 0.25) 


2.0 - 2.6 


ko 


<£, fc . > = (0.33 - 0.44) xlO" 2 


2.0 - 3.5 


(£ 2 , fc .) = (0.352 - 0.358) 


1.8 - 2.3 


fo 


itlfn) = (0.331 - 0.339) 


1.8 - 2.2 


(Clfn) = (0.204 - 0.220) 


1.8 - 2.2 



\Ss,a / 



0.369 (0.245), 
= 0.004 (0.355), 
= 0.335 (0.212), 



1(2) 
A"* 



(f 1(2) 

(Qf ) = 0-196 (0.088), 



0.203 (0.093), 
= 0.018 (0.207), 



(23) 



It is noted that these moments are different from Ref . 0] , 
especially for the moments of Kq and fo. By using the 
same input parameters, we can obtain consistent mo- 
ments for ao Then, the differences are mainly caused 
by the different treatment of the mass terms (the s-quark 



By using the masses and decay constants for the 
scalar mesons, we further calculate the moments of these 
mesons from the sum rules (0 UHl HU H2J • Considering 
the conservation of charge parity and isospin symmetry, 
the odd moments for ao and fo meson twist-3 distribution 
amplitudes should vanish. While taking SU(3) symme- 
try breaking effects into account, the odd moments for 
Kq meson are nonzero. Therefore, the moments we are 
aimed to derive are the second and fourth moments for oq 



-«..*;} 



-(&«;> 
-(&«,;> 



FIG. 5. Moments of scalar meson Kq versus M 2 within their 
Borel windows. 
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-<&*> 



-(ft/.) 

-<a,/.> 



FIG. 6. Moments of scalar meson /o versus M within their 
Borel windows. 



mass terms) involved in the hard part calculation 0. Sim- 
ilar to the kaonic case Q , it is found that those m s -terms 
do provide sizable contributions. Thus, they should be 
treated consistently with those of the higher dimensional 



matrix elements. By taking all the mass terms consis- 
tently into consideration, a more reliable masses and de- 
cay constants, and hence more accurate moments can be 
obtained. 



A discussion on the scalar's three particle 
twist-3 LCDAs 



For each scalar meson S, there are three twist-3 distri- 
bution amplitudes 4> s s , 4>g and (j> 3 s- To be useful refer- 
ence, we make a discussion on the three particle twist-3 
LCDA, which, similar to the pseudoscalar and the vector 
cases [lalMliJ], is defined as 



\q2(x)a^gG aP {-vx) qi (-x)\S{q)) = if 3S [q a (q^ - q u S^) - (a O /?)] / Va t e ^(-«i+«+«a 3 )^ 3S ( ai ) ( 24 ) 



where T>oti — da\da2da 3 b{a.\ + a 2 
Ri 



1 ha D 1 f3K* 
• '>2 



m a fa 



m K * f K * 



a 3 

Rs = 



1). Defining 

1 /g/o 
l h ffo ' 



by using the recurrence relations for the moments among 
the twist-3 LCDAs that are derived by using the equa- 
tions of motion and the conformal symmetry (43l - l45{ 0, 
we obtain 



(£cr,ao / 



g<0 + f «l-g^««3)} aQ 



2 ) 
s,ao / 



(t,s,K * ) 



4i?i 



— R 3 

5 

4i? 3 
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^ 3 ((a 3 ))/o 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 



where (((a 2 - on + va 3 ) n )) = J T>a i 4> 3 s{ct i )(a 2 - ai + 
va 3 ) n defines the moments of the three-particle twist-3 
LCDAs. Using the central values for the two particle 
twist-3 LCDA moments, we obtain 

«a 3 »a ~ 1-29 and f 3ao ~ 4.82 x 10- 3 GeV 2 , (31) 

((as))*. ~ 2.08 and f 3K * ~ 2.76 x 10- 3 GeV 2 ,(32) 

«"3)>/ ~ 7.50 and f 3ao ~ 2.56 x lO^GeV 2 . (33) 



1 Since the current s-quark mass is at the order of Aq^q, it is 
reasonable to assum that those terms will have sizable effects. 

2 As an estimation, we neglect the meson mass effects in deriving 
those relations. 



These equations show that all Ri are less than 1%, so 
in usual calculations, one can safely neglect the scalar 
meson's three-particle twist-3 LCDA in comparison to 
the two-particle twist-3 LCDAs. This is different from 
the pion and the kaon, whose similar ratios R^ {Rk) 
~ 5% [45j : then for the pionic or kaonic processes, one 
may need to take (f> 37r or <f> 3 K into consideration. 



D. Scalar meson distribution amplitudes 

The scalar mesons' twist-3 LCDAs can be expanded 
into a series of Gegenbauer polynomials jZJ, |42| 

00 

^( U ,^) = l+^a m ( At )Ci/ 2 (2 U -l), (34) 

m— 1 

00 

n) = 6uu[l + b m (v)CV 2 (2u - 1)]. (35) 

m=l 

where Cm 2,1 ^ 2 (2u — 1) are Gegenbauer polynomials. The 
Gegenbauer moments a m , b m can be related to moments 
and defined in Eqs.jSl H- By using the or- 

thogonality of Gegenbauer polynomials 

1 duC]l\2u - l)Cl' 2 {2u - 1) = -^—S mn , (36) 

In + 1 

1 

duu(l - u)Cl /2 {2u - l)CH 2 {2u - 1) = 



(n + 2)(n + l) 



4(2n + 3) 

we obtain the usual Gegenbauer moments: 



ai = 3<6), a 2 = - (3(6) - 1) , 



(37) 



a 4 = - (35(&) -30(6) +3), 



(38) 



&4 = ^ (21(6) - 14(6) + 1) 
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(39) 



TABLE V. Gegenbauer moments for <f> s s at two energy scales 1 GeV and 2.4 GeV, where S stands for the scalar meson oo, kg 
or fo respectively. 



mesons 


H = lGeV 


H = 2.4GeV 




ai 


a 2 




ai 


ai 




ao 





0.302 ~ 0.323 


0.491 ~ 0.944 





0.225 ~ 0.241 


0.307 ~ 0.591 


kg 


0.0109 ~ 0.0143 


0.163 ~ 0.211 




0.0095 ~ 0.0125 


0.122 ~ 0.157 




fo 





-0.022 ~ 0.049 


0.140 ~ 0.892 





-0.016 ~ 0.036 


0.088 ~ 0.56 



TABLE VI. Gegenbauer moments for (j>g at two energy scales 1 GeV and 2.4 GeV, where S stands for the scalar meson ao, kg 
or fo respectively. 



meson 


H = 1 GeV 


At = 2.4 GeV 


bi 


b 2 


b 4 


bi 


62 




ao 





0.011 ~ 0.013 


0.050 ~ 0.073 





0.008 ~ 0.009 


0.031 ~ 0.045 


kg 


0.0216 ~ 0.0468 


0.019 ~ 0.028 




0.017 ~ 0.037 


0.013 ~ 0.020 




fo 





-0.0291 ~ -0.0014 


0.030 ~ 0.106 





-0.0206 0.0010 


0.018 ~ 0.065 



Using the moments and derived in the last 
subsection, one can obtain the Gegenbauer moments a m 
and b m by using the relations (|38I39[) . For convenience, 
we present the first two non-zero ones for each meson at 
two different energy scales in Tables [V] and IVI1 where 
the Gegenbauer moments at fi = 2.4 GeV are obtained 
by using the evolution equations. 

The LCD As 4> s s {u) an d <f>g(u) at fi = 1 GeV are pre- 
sented in Figs. (I7I8|) . As a comparison we also present the 
results of Kef. 2 in Figs. (011). 

E. Properties of the B — > S transition form factors 

In doing the numerical calculation we adopt 0, H3, 13 : 
m Bo = 5.279 GeV, m Bs = 5.368 GeV, f Bo = (0.19±0.02) 
GeV, f Be = (0.23 ±0.02) GeV, m b = (4.68 ±0.03) GeV. 
As for the energy scale of the B — > S transitions, we 

adopt /ib — J m Bs — m 1 ~ 2.4 GeV. To calculate the 

form factors one has to evolve the parameters to the scale 
/ib, which can be derived by using the renormalization 
group equations (IA1IA2I) . We take the threshold param- 




FIG. 7. Twist-3 LCDAs 4>s( u ) f° r scalar mesons at the scale 
fi — 1 GeV, where the solid, the dashed and the dotted lines 
denote the LCDAs of ao, Kq and fo respectively. The left 
diagram is for our LCDAs and the right one is for Ref.[2|. 



eter S* B ° = S Bs = (33 ± 1) GeV 2 30 . The Borel window 
is determined by requiring both the contributions from 
higher excited resonances and continuum states be less 
than 30%, which results in 10GeV 2 < M 2 < 15GeV 2 . 
As shown in Fig. (J9j) , the sum rules for the form factors 
at the large recoil region vary slightly within such Borel 
window. 

We present the B — > S form factors at the large recoil 
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TABLE VII. Numerical results for the B —¥ S transition form 
factors at the large recoil point q 2 = 0, where S stands for 
the scalar mesons a + , fo, Kg + and Kq, respectively. 




FIG. 8. Twist-3 LCDAs (/>'s( u ) f° r scalar mesons at the scale 
fi — 1 GeV, where the solid, the dashed and the dotted lines 
denote the LCDAs of ao, Kq and fo respectively. The left 
diagram is for our LCDAs and the right one is for Ref.0j. 



M"(GeV ) 

FIG. 9. Prediction of f+(q 2 = 0) within the LCSR approach 
with chiral currents versus the Borel parameter M. The solid, 
the dashed, the dotted and the dash-dot lines are for Bo — > 
a+(1450), So ->■ i?o(1430), B° ->■ #*+(1430) and B° -> 
/o(1500) respectively. The input parameters are taken as their 
central values. 



processes 


Methods 


£ 

f+ 


£ 

I- 


£ 

JT 


B° - 


■¥ a+(1450) 


This work 


n /i zi+0.06 




0- 4 3I .o5 






LUan 140] 


n co 


f\AA 
— U.44 


U.DD 






LGSK 130] 


0.53 


—0.53 








PQLjU 1471 


U.DO 




c\ no 






PIT I A Ql 


U.zo 






s° - 


+ /o(1500) 


This work 


A QO+0.04 
O.J8_ .04 


n 0,1 +0.04 


n , n +0.04 

0.40to.o4 






LljOll 140] 


U.4o 


— U.O f 


U.OD 






t pcb ronl 
LUbK IoU) 


n A 1 
U.41 


n /i i 
— U.41 


u.oy 






p^OD |47] 


U.ou 




n on 

U.OZ 




■¥ #*+(1430) 


This work 


n QQ+0 04 


n O r+0.05 
U - /0 -0.05 


n +0.04 
0' 41 -0.04 






T nnn [AH] 

LUbK |4oJ 


U.42 


— U.34 


U.OZ 






t noD ronl 
LUbK IoUJ 


U.44 


— U.44 








an 149] 


0.24 










pQCD [47] 


0.56 




0.72 




+ iv *(1430) 


This work 


U.4O_ 05 


n oa +0.06 


4fi+ 06 






LCSR [46] 


0.49 


-0.41 


0.60 






LCSR [30] 


0.49 


-0.49 


0.69 






SR [50J 


0.31 


-0.31 


-0.26 






pQCD [47] 


0.60 




0.78 






CLF [48J 


0.26 










LFQM [51] 


-0.26 


0.21 


-0.34 



point q 2 — in Table IVII1 where the theoretical uncer- 
tainties are around 20%, which are caused by varying the 
Borel parameter M, the threshold values Sq° and Sq* , the 
b quark mass, the decay constants of the scalar mesons, 
and the Gegenbauer moments for the twist-3 LCDAs of 
the scalar mesons, within their reasonable regions. As a 
comparison, we also put several typical estimations from 
the SR, the LCSR and the pQCD approaches USEUm, 




q 2 (GeV 2 ) 

FIG. 10. Form factors f+(q 2 ) within the LCSR approach 
with chiral currents. The solid, the dashed, the dotted and 
the dash-dot lines are for Bo -> (4(1450), B -> #3(1430), 
B° -¥ K* + (U30) and B° -¥ / (1500) respectively. Here the 
Borel parameter M 2 — 12 GeV 2 within the LCSR approach. 




FIG. 11. Differential decay width for the B —> Slvi decays. 
Here the massless lepton I — e or fi. Here S stands for 
a + (1450) or # + *(1430). 



in Table IVlTl The form factors f+(q 2 ) for various scalar 
mesons a^(1450), #£(1430), ^+(1430) and / (1500) 
rise slightly with the increment of q 2 , as can be shown 
explicitly in Fig. ([TU)) . 



F. Two semi-leptonic channels B — > Slui and 
B -+ Slut 

Next, as an application of the B — > S transition form 
factors, we discuss the properties of the two semi-leptonic 
channels, B — > SWi and B —> Slui, whose differential 
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FIG. 12. Differential decay width for the rare B — > Sll de- 
cays. Here the massless lepton I — e or fi. Here S stands for 
7? *(1430), /o(1500). 



decay widths can be expressed as 



dT G 2 F \V u b\ 2 q 2 -m 2 

— [B {s) -+ blvt) - mn3m3 {q2)2 



(q 2 — mf) 2 / (m B — m s — q 2 ) 



2\2 



(m 2 + 2q 2 ) 



(q 2 - [m B - m s ) 2 )(q 2 - (m B + m S ) 2 ) fl{q 2 ) + 



3mj(m B -m 2 s y [ f + tf) + — 
1 mi, 



,(40) 



and 



dT 
dq 



— (B {s) ^Sll) = 



G 2 F \V tb V ts \ 2 m B a 2 em 



4n 

s 



1/2 



1536tt 5 

1 , 2ri \ 3 / 2 , V2 x 



(41) 



where m/ denotes the mass of a final state lepton, and 



2 / 2 2 / 2 2 / 

s = q /m B , n = m l /m B , r s = m s /r 
<ps = (l-r s ) 2 -2s(l + r s )+s 2 , 



B- 



(42) 
(43) 



a s 



r eff f ( n 2\ g 7/r(g 2 ) 



|C 10 /+(<Z 2 )|!44) 



5s-6|C 10 | 2 {[2(l + r s )- S ]|/ + ( (7 2 )| 2 + 

(1 - rs) 2Re[f + (q 2 )r_(q 2 )} + s |/-(g 2 )| 2 }. (45) 

Except for the B — > S transition form factors, we adopt 
the same input parameters as those of Ref. [30| ■ Our re- 
sults for the curves for the differential decay rates are 
presented in Figs.lfTTl [T2"|) . which have a different aris- 
ing trends versus the increment of q 2 in comparison to 
that of Ref. [3(| . In Ref. [30[ , the authors also adopt the 
chiral currents in calculating the B — > S transition form 
factors. However their chiral currents are different from 
ours, where only the twist-2 scalar LCDAs are remained 
and the twist-3 scalar LCDAs make no contributions. 
Those two different QCD sum rules can be taken as a 
cross check of the chiral LCSR approach and be con- 
firmed by the future experiments. 



IV. SUMMARY 

The masses, the decay constants and the twist-3 
LCDAs of the scalar mesons do, Kq and /o have been 
investigated in this work by using the QCD sum rules 
together with the QCD background field theory. Our 
present estimations are improved by a consistent treat- 
ment of the mass effects and a RG improved treatment 
of the input parameters. It is found that the m s -terms 
provide sizable contributions, which should be treated 
consistently with those of higher dimensional matrix ele- 
ments. While, by taking all the mass terms consistently 
into consideration, a more reliable masses and decay con- 
stants, and hence more accurate moments can be ob- 
tained. 

As for do meson, the second and forth moments of the 
twist-3 LCDAs ^ are 

(& = 0-369 (0.245) and (£ 2 ^) = 0.203 (0.093), 

whose uncertainties are about 10%; as for the Kq meson, 
the first and second moments of the twist-3 LCDAs 6 S ^Z 
are 



= 0.004 (0.355) and = 0.018 (0.207), 

whose uncertainties are about 10% — 15%; as for the case 
of /o meson, the second and forth moments of the twist-3 



1(2) 



= 0.335 (0.212) and (£<j>) 



= 0.196 (0.088), 



whose uncertainties are within 15%. 

As an application of these twist-3 LCDAs, we have 
studied the B — ^ S transition form factors and their 
corresponding decay rates. For the purpose, the chiral 
currents are adopted in our LCSR calculation of B — >• S 
transition form factors, in which the twist-3 LCDAs make 
dominant contributions. Our results for the transition 
form factors at the large recoil point q 2 ~ are con- 
sistent with those obtained in the literature; while, the 
arising trends for the form factors versus q 2 , and hence 
the differential decay widths for the B —> S semileptonic 
decays, are somewhat different. This can be checked by 
a future experimental data. 

Our results for the twist-3 LCDAs are useful for 
pQCD calculation. Basing on the Gegenbauer moments 
of the twist-3 LCDAs, if further taking a proper trans- 
verse momentum dependence, e.g. the BHL-transverse 
momentum dependence (52j, we can construct a better 
behaved scalar wave functions in which the end-point 
singularity can be effectively suppressed. For example, 
following the similar idea in constructing the pseudo- 
scalar meson's twist-3 WF model, it is natural to assume 
that one can obtain reasonable twist-3 contributions to 
the scalar meson involved high energy processes. 
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Appendix A: Details for deriving the twist-3 LCDAs 

The LCDAs of the scalar mesons are normalized as 
J du(f) s s (u, fx) = J Q du4>g(u, /i) ~ 1. The initial energy 
scale for the bound state and the non-perturbative matrix 
elements are set as the values of the Borel parameter M. 
Their values at any other scale /i can be obtained from 
the RG equations, for example, the decay constant, the 
quark mass and the condensates are related by (3(J [4(| 



fs(M)=f s (p) 



m q (M) = m q (n) 



4/6 



a s (M)J 



-4/6 



(qq) m = (qq) f 



4/6 



{g s qaGq) M = (g s q<rGq) 



a s {M)J ' 

, \ x -2/36 



a s (M) 



>M 



(Al) 



where b = (33 — 2n/)/3 with nf the active flavor number. 
And the RG equations for Gegenbauer moments are [53| 



(a n (fi)) = (a n (fj, )) 



o 



M X 



-illh 

where the one-loop anomalous dimensions are 

n+l 



>£=Cf 1- 



8 4 V- 

(n+l)(n + 2) + j 

3 — A 



n+2 



• t -°'( 1+4 £j 



J'=2 



(A2) 

(A3) 
(A4) 



with C F = 4/3. 

The moments of the LCDAs for the scalar mesons are 
calculated under the background field approach. For the 
purpose, we adopt the following two correlation functions 



i J dW^OllWsXi* Drq 2 (x),q 2 (0) qi (0)}\0) = -(z-q) n lt \q 2 ), (A5) 
xe^'WTfaWa^iz- D) n+1 q 2 (x), &(0)?i(0)}|0) = - q v z^(z ■ q) n I^\q 2 ). (A6) 



On the one hand, the correlation functions (|A5IA6p 
can be evaluated using operator product expansion in the 
deep Euclidean region (-q 2 > 0). With the help of the 



higher-dimensional mass terms reserved in the expression 
of the propagators ([7|) and |[8} , the operator product ex- 
pansion can be expressed as 



4V) ocfl ^ dx(2x-l 



QCD - 4?r 2 

+ {n + 2)[q 2 x(x-l) + 



) n ^[q 2 x(x- 1) +mim 2 ] ML -In 



q 2 x(x — 1) + m 



L 12\ 





1 dxx(l 
'0 


^ n(n-l) f 




qiqi) 


2m 2 + (1 


2(m| 



3 q 2 x(x — 1) + m\ 2 
2 ~ n ~a~ 2 



n+l 
5(21 — 1) H 
l (l-x) 2 



2mim2 — mf 2 



,2 12 
i 12j 



i-2 



+ - 



mfmf 



{g s qiaTGqi) 



6(m,2 — q 2 ) 2 [vn\ — g 2 ) 3 

[18(71 — l)?7i2 + n(8n — 5)mi] [(An — 3)mi + 3m2\m\ 



36(m§ - <? 2 ) 



6(m2 — g 2 )' 



13 



and 



2m\m 



3K - <7 2 ) 4 
-9s (9i ?i) (9292) 



5s(?i9i) 2 



2n 2 + 7n - 12 



2ml 



7 2 + 2mi?7i2 



9 g 2 (m 2 — q 2 ){m\ - q 2 



Sl{m 2 -q 2 ) 2 27(m 2 -g 2 ) 3 

+ (-l) n (t?i <-» g 2 ,mi O m 2 ) 



9 2 s(liqi)(q2q2 



4 to 2 + g 2 
9 <7 2 (m§ — q 2 ) 2 



3( " + 9 1} ^(2* - l)Vs(* - 1) + m 2 , 
4tt 2 Jo 

1 " + ^-<a s G 2 ) / da:(2:r - l)"x(l - &) 



3 j q 2 x(x — 1) + m 



24tt 



. (ff 3 /G3) (- + l)-i- 2 f 1 



19191/ 



24tt 2 
-mi(n + 1) 
6 



/ dz(2a-l) n 
Jo 



2 ^ 

q 2 x(x — 1) + m 2 2 — 2mim2 
[g 2 x(a; — 1) + m 2 2 ] 2 

,t 2 (1 -x) 2 



[g 2 x(a: — 1) + rrc 2 
mf(2n + 1) 2m\m 2 , 



2 13 



+ (g s qiaTGqi){n + 1) 



+ff 2 (9i9i)V+l)[ 



2 2 — 5 2 ( m 2 ~~ 9 2 ) 2 ( m 2 — 9 2 ) 3 
3n + l)mi + 6m2 2rti\m\ 



36(m 2 - 9 2 ) 2 9(m 2 -g 2 ) 3 



-2n + 5 



4 m, 



81(m 2 -g 2 ) 2 81(m 2 - g 2 ) 3 



+ (-l)"(<?i <-> g 2 ,mi O m 2 ) L 



(A7) 



(A8) 



where 



l 12 



to 2 x + m 2 (l — x). 



I 

can be derived by inserting a complete set of quantum 



On the other hand, the correlation functions (|A5IA6p states S|n)(n| in the physical region: 
I 



Im/ S "(g 2 ) w = ~7rS(q 2 - m 2 )m 2 f 2 (C) + ^ / dx(2x - l) n \ (n + 3)q 2 x(x - 1) + (n + 2)m\ 2 + mi m 2 \d(q 2 - s s ) 



n +!„ 2 52/cnv , 3(n + l) Z" 1 



4?r 



dx(2x - l) n ^ q 2 a;(x - 1) + m 2 2 ^(q 2 - S<r ) 



In deriving the above equations, we have implicitly 
adopted the quark-hadron duality. These two expres- 
sions of correlators (|A5|A6j) can be matched through the 
dispersion relation 



ds 



lml(s) had 2 

— = * [q )qcd- 

s — q z 



(A9) 



One can apply Borel transformation to both sides to sup- 
press the unknown higher dimensional condensates and 
the continuum contributions as much as possible. For an 
arbitrary q 2 polynomials, the Borel transformation im- 
plies [Hi 



B M 2 = lim 



-q 2 /n=M 2 



r 



where M is the Borel parameter. 



Applying the Borel transformation to the dispersion 
relation, it results in 



1 J dse- s / M2 lmI(s) had = B M 2l(q 2 ) QCD , 

which finally leads to the sum rules P [TT1 ITUl \12~\i 

adopted in the body of the text. 
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